Lagrangian Formulation of Stochastic Inflation: A Recursive Approach 
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We present a new, recursive approach to stocliastic inflation which is self-consistent and solves 
multiple problems which plagued a certain number of previous studies, in particular in realistic 
contexts where the background spacetime is taken to be dynamical, where there is more than 
one field present, especially with a mass hierarchy, or where the role played by back-reaction is 
suspected to be important. We first review the formalism of stochastic inflation as it is usually 
heuristically presented, that is, deriving the Langevin equations from the field equations of motion, 
and summarize previous results on the subject. We demonstrate where inconsistent approximations 
to the Langevin equations are commonly made, and show how these can be avoided. This set up 
shares many similarities with quantum Brownian motion and out-of-equilibrium statistical quantum 
dynamics. We hence review how path integral techniques can be applied to the stochastic inflationary 
context. We show this forinalism to be consistent with the standard approach, develop a natural 
perturbative expansion, and use it to calculate the one-loop corrected Langevin equations. 
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verse, it becomes more and more critical to develop solid 
and self-consistent methods of calculation for inflation- 
ary predictions. For example, in the context of multi- 
field inflation, it is complicated to disentangle the grav- 
itational and matter degrees of freedom when describ- 
ing fluctuations produced in the scalar fields using tradi- 
tional methods. Typically, approximations are made to 
make the problem tractable which ignore back-reaction, 
that is, the effects of these very fluctuations on the back- 
ground spacetime and fields trajectory. Restoring or even 
assessing the importance of these neglected effects then 
becomes extremely non-trivial. 

For instance, in [5[, we presented an analytical study 
of mode coupling in the context of hybrid inflation. More 
precisely, we have studied the effects of short wavelength 
fluctuations both on the background inflaton and on long 
wavelength modes of the inflaton and waterfall fields. 
However, we found that our results were strongly depen- 
dent on initial solution chosen, as well as on the choice of 
background. Because of this sensitivity, a more thorough 
study of the background was necessary, a goal which can 
be achieved via a study of stochastic inflation [6l-ll4i| . In- 
deed, in [15|, it was shown that stochastic effects can 
significantly alter the inflationary background dynamics 
in the context of hybrid inflation. 



I. INTRODUCTION 

Since the publication of the seminal paper by Guth 
[l|, where the idea of inflation was flrst proposed, infla- 
tionary cosmology has become one of the cornerstones 
of modern cosmology. Indeed, not only did it propose 
an elegant solution to both the homogeneity and flat- 
ness problems, but being the flrst compelling cosmolog- 
ical model allowing predictions based on causal physics 
about the structure of the Universe on large scales, this 
theory has dramatically changed the face of cosmology. 
The predictive power of inflation [2[ together with the ob- 
servations of the Cosmic Microwave Background (CMB) 
anisotropics with the COBE and the WMAP H exper- 
iments, and very recently the Planck satellite [J|, have 
truly embodied the last step in this transformation. Fur- 
thermore, by opening new avenues connecting theoretical 
physics to experiment, it has allowed substantial progress 
in ways to approach fundamental high energy physics. In 
fact, through the prediction of the shape and amplitude 
of the spectrum of primordial fluctuations of the CMB, 
inflation has provided us with the theoretical tools to 
analyze what can be thought of as a direct window into 
the Planck scale. In this sense, CMB experiments give 
us a direct access to the quantum world at energy scales 
far beyond what is accessible in Earth-based high energy 
particle physics experiments. 

With the ever-increasing precision of the experiments 
probing this window into the very beginning of the Uni- 
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The stochastic inflation formalism is very powerful 
when it comes to addressing such issues, because it al- 
lows for a continual renormalisation of the background 
dynamics. Its strength lies in its ability to separate the 
dynamics of coarse-grained long, classical, super-Hubble 
wavelengths from short, quantum fluctuation-dominated 
sub-Hubble wavelengths in a way that constantly cor- 
rects the background dynamics as modes of fluctuations 
are stretched from the quantum regime into the coarse- 
grained effective theory. The theory then describes the 
effective classical dynamics of the large-scale degrees of 



freedom, in the presence of a "bath" where all the quan- 
tum fluctuations are collected in a classical noise term. 
Rather than being simply a computational trick, this fea- 
ture is claimed to be fundamental to the system under 
study [7[ . One of the major advantages of this formalism 
is then that, the super-Hubble theory being classical, one 
can solve it non-perturbatively, or to arbitrary order in 
perturbation theory [16|, whereas it is unlikely one will 
ever be able to solve exactly the Heisenberg field equa- 
tions for four dimensional interacting quantum fields. 

The resulting theory shares a lot of similarities with 
quantum Brownian motion and out-of-equilibrium quan- 
tum statistical mechanics. Indeed, the classical system 
consisting of the super-Hubble modes is in contact with 
a bath of quantum fluctuations, modelled as a random 
noise source which must be averaged over time. The sta- 
tistical properties of the noise are then all that is needed 
to derive the probability distribution of the system's state 
over many realizations. 

This link with out-of-equilibrium mechanics suggests 
that one can derive the formalism to study inflationary 
dynamics from the machinery usually applied to quan- 
tum Brownian motion, such as the influence functional 
formalism [17| . It arises in the context of the Schwinger- 
Kcldysh closed time path (CPT) formalism, also known 
as the in-in formalism, which has become widely applied 
in cosmology since [la Um i after the pioneering work of 
[20I [2]| in the 1980's. However, in the case of the in- 
fluence functional formalism, the strategy is to split the 
degrees of freedom of the full quantum flelds into mo- 
mentum space through a window function, and perform 
the path integral over the small scale fluctuations. 

In the present paper, we will flrst derive the stochas- 
tic formalism in Section II starting from the equations of 
motion of the fields, as was done in the early work of [3- 
[12. 22 - 28,j . However, being very careful in the derivation, 
we obtain terms which are not usually kept and we see 
precisely at what level the standard approach is making 
crucial, and often inconsistent, approximations. With 
this understanding, we present a new method to con- 
sistently treat the problem of back-reaction and mode- 
coupling in the context of stochastic inflation. 

The key issue is that, during inflation, the propagator 
of the bath degrees of freedom evolve in a background 
which has to be flxed in order to compute the amplitude 
of the noise sourcing the bath, coarse-grained degrees of 
freedom. However, this effect shifts, or renormalizes, pre- 
cisely this same background, making the propagator shift 
again. That is, as we will see, the bath degrees of free- 
dom depend on the physical, renormalized background, 
not on a fixed free value as has been assumed in past 
studies of stochastic infiation. This makes the problem 
of integrating out the bath degrees of freedom much more 
difficult. Indeed, one must now make sure to evaluate at 
the same order both the coarse-grained flelds and the 
quantum flelds. We propose a new, recursive method to 
achieve this precise goal at the end of Section II. 

In Section HI, we explore a more rigorous approach in- 



spired by condensed matter physics and prove that both 
derivations yield the same result for the Langevin equa- 
tions at linear order. As we will see, one of the major 
advantages of the so-call influence functional formalism 
is that it makes it very natural to include higher loop 
corrections to mode-couplings and back-reaction effects. 
With this in mind, we develop a natural perturbative 
expansion and use it to calculate the full one-loop cor- 
rected Langevin equations in section IV, which is, to the 
best of our knowledge, a new calculation in the context 
of inflationary (with a time-dependent window function) 
spacetimes (see, however j29ll3G|). 

With this machinery developed, we then have all the 
necessary tools to apply our new recursive method to the 
speciflc example of two-fleld hybrid inflation potentials. 
We shall do such thing in a second paper |3l| . 



II. STOCHASTIC FORMALISM: A HEURISTIC 
DERIVATION 



We would like to describe the evolution of scalar fields 
minimally coupled to gravity in de Sitter space. This 
task is very difficult to perform accurately using semi- 
classical methods, since one would expect effects such 
as self-coupling and interactions to become more and 
more important as time passes and the phase space of 
modes which have frozen and which are growing outside 
of the Hubble horizon becomes correspondingly larger 
and larger. To tackle this problem we take advantage of 
the similarity of the current system with quantum Brow- 
nian motion. That is, we wish to use the techniques of 
the well-established Schwinger-Keldysh theory for out-of- 
equilibrium quantum field theory to derive the stochastic 
nature of the classical limit, in the form of a Langevin-like 
equation with a non-linear system-bath coupling. The 
techniques used in this case can be very technical and can 
even obscure the physical meaning of the approximations 
made. Therefore, we first present a heuristic derivation 
from the equations of motion, before proceeding to the 
more formal quantum field theoretic path integral ap- 
proach. 



A. Definitions and formalism 

We consider scalar fields minimally coupled to gravity 
and evolving in a spatially flat background Fricdmann- 
Lemaitre-Robertson- Walker space with the metric: 



ds^ = -dt^ 



\t)d2 



^(r)(-rfr2 + rff2) ; (1) 



where the second equality links the conformal time r with 
the cosnuc time t hy dr = dt/a{t). Derivatives of a(t) 
give the Hubble parameter, H , and the deceleration pa- 



ramctcr q{t): 



"'"'-ly -ik' 



, ^ aa H 



(2) 

Here, we have related H to the matter energy density 
represented by p through the first Friedmann equation. 

For inflation to occur, we require positive expansion 
{H{t) > 0) with -1 < q{t) < (where the lower bound 
comes from the weak energy condition) . When the decel- 
eration q saturates the lower bound, H is constant and 
we obtain a locally de Sitter geometry, in which a(t) is 
exactly exponential. In what follows, we shall not fix the 
background geometry to be exact de Sitter space, and 
consider quasi-de Sitter space (in the sense that we allow 
for H y^ 0). Additionally, for the sake of clarity of the 
notation, we will consider an unperturbed homogeneous 
and isotropic space-time in the current paper, and post- 
pone the inclusion of metric perturbations to the second 
paper of this series, where we implement our method in 
a concrete example. 

We are interested in studying effects due to the pres- 
ence of multiple fields during inflation, and will therefore 
introduce a formalism treating two scalar fields $ and ^. 
Generalizing it to a larger number of matter fields should 
be straightforward. The matter sector of the action we 
are considering is described by: 



with 

/:,„($, *) = 



dt I d-'^x 



-lid-i>r~lid^r~vi^,^) 



(3) 



(4) 



It will be useful to split the potential into mass terms 
and interaction terms: 

Vi<^, *) = im|$2 + im^*2 ^ Vp,,ti^, *) . (5) 



B. Coarse-grained and quantum fields 



In the spirit of the early work on the subject [7H12I.I22 
l28l |. we derive here the stochastic formalism using the 
fields equations of motion as a starting point. The un- 
derlying idea is to study them in Fourier space and split 
the modes into sub-Hubble and super-Hubble. We will 
be interested in studying the effective classical dynamics 
of the large-scale fields by treating the sub-Hubble modes 
collectively (through a coarse-graining process which will 
be defined shortly) as a classical noise source. 

Our strategy is the following: we wish to consider 
the full quantum fields ^ and $, and split them into 
a large (in both amplitude and scale) "homogeneous" 
piece and an inhomogeneous piece representing small, 
high wavenumber, quantum fluctuations: 



where, by "homogeneous", we mean the full flelds <&, ^, 
coarse-grained over some large enough volume. Strictly 
speaking, it is an oversimplification to think about the 
coarse-grained fields as classical and it would be more 
accurate to denote (p and x by quantum mean fields (af- 
ter all, from what has been said, both the coarse-grained 
and inhomogeneous fields are quantum by nature) . How- 
ever, as we will see later, it turns out they can actually 
be identified with the classical mean fields, while </>> and 
ipy can be identified with the quantum fluctuation fields. 
Also, see [l2, 120, l32| - |43 | and references therein for a dis- 
cussion of decoherence and the precise condition required 
for quantum fields to undergo classicalization. 

The question of defining the coarse-graining volume in 
a precise manner might also seem non-trivial, and we will 
quickly discuss it in a moment. For now, let us simply 
note that choosing an invariant scale seems like a sensible 
thing to do (in the same way that it is crucial to choose 
an invariant ultraviolet regularization scheme in order to 
regularize, for example, scalar fields UV divergences in 
the context of general relativity) . 

To give us a hint of what this invariant scale should be 
chosen to be, consider a given fixed comoving wavenum- 
ber k = |k| of a free, massless field sitting in de Sit- 
ter space. For this mode, significant particle production 
starts to occur when k < Ha{t). Moreover, for free, mass- 
less scalars in exact locally de Sitter space, the range of 
physical wavelengths from zero to a constant H^^ will 
have a constant dynamical impact since they lie in an in- 
variant range [l6| (provided we choose an invariant reg- 
ularization prescription). These facts all point toward 
using the de Sitter length as a cutoff scale to separate 
the coarse-grained fields (p and x from the quantum fields 
(/)>, ')/'>. However, since we won't be considering exact 
de Sitter space in the following, and won't be restricting 
ourselves to massless fields, we expect that the above will 
only be approximately true in the case under study here. 
Regardless, this still suggests a natural choice of physical 
separation for the coarse-grained fields ((y9, x) 3'iid ultra- 
violet fields ((/)>, '0>), so that, in terms of the Fourier 
modes k of the fields. 



> , ■0> correspond to 
</? , X correspond to 



k > H{t)a{t) , (7) 

H(t)a{t)>k>0, (8) 



C. Effective classical coarse-grained dynamics 

The equations of motion satisfied by the full fields are: 



n$ + m|$ = -Vrt,* (*,*), 



-n* + mi* = -K 



q,^ - -•'pcrt,*($,*) 



(9) 
(10) 



$ = V' + 0>, ^^x + V"; 



(6) 



where —0 = dtt + 3Hdt — V^/a^ is the d'Alembertian 
and the commas denote a partial derivative(s) with re- 
spect to the field(s) following it. Inserting ^ into © 



and assuming that the quantum fields are a smah per- 
turbation to the coarse grained field, i.e. Lp ^ (/>> and 
X ^ V'>j we write: 



-n 



-□(/? + m%ip + Vp„,,<j,((^, x) 
<P> + ml(f>y + V^'^ {if, x)(f>> + ^$* {(f, x)'0; 



= -yS,*',T,(<p,x)(?!>>'0; 



,***>. 



-hv:^'^^i^,x)<Pi 



1 T rpcrt 

2 ''^,*** 



(¥',x)V'| + 



(11) 



where the ellipses stand for terms of third order or greater 
in the quantum fields. We obtain a similar equation for 
\E'. The second line of the above equation, that is, the 
one in brackets, corresponds to the linearized equation 
satisfied by the small quantum perturbations 0>, in the 
presence of the mean field (p background. Subtracting it 
from the full equation, one might think that one would 
find the equation for the coarse-grained field. 

However, things are slightly more complicated since 
the cutoff in Fourier space defining the comoving 
wavenumber range over which the quantum field extends 
is time-dependent. To see this, we start by expanding 
0> and ipy into creation and annihilation operators on a 
time-dependent background: 






V'k^ke 



-ik-x 



t ^ik-x 






^t 



f,t 



and &k, ^k obey the standard 



where the operators Sk, 
commutation relations. 

Here Wi/(fc,i) is a window function that filters only 
the high-fc modes and projects out the long-wavelength 
modes with k < a{t)H{t). The simplest choice of win- 
dow function is the theta function 9 [k/{eaH) — 1], where 
e < 1 is some constant small enough to ensure we include 
all the modes with a genuinely oscillatory quantum na- 
ture in the two quantum fluctuation flelds. However, as 
was shown in |43| . a sharp cutoff in k-space might be 
somewhat pathological, even if one could argue that such 
a sharp cutoff is only the mathematical limit of a smooth 
filtering procedure. Indeed, it has been shown [ij] that 
it does not satisfy the same fundamental properties as 
the asymptotic limit of any of these smooth limiting pro- 
cesses, and therefore has the physical effect of biasing 
the correlations of the coarse-grained fields on the largest 
scales to spuriously high values. 

Most importantly, when the quantum fields are treated 
as free, massless fields in de Sitter space, it is the only 
window function that will cause them to appear as white 
noise to the coarse-grained fields, which means the latter 
will behave as a so-called Markovian process [131 • This 
means that at every step of its evolution, the system 
made of the coarse-grained fields will experience noise 
which has no memory of previous history. On the con- 
trary, considering a well-behaved smooth window func- 



tion gives rise to colored noise, and consequently non- 
Markovian processes. Even though their treatment is 
more challenging, they are more physically motivated 
since they correspond to a finite volume cutoff between 
quantum and coarse-grained fields in configuration space, 
rather than a complicated infinite volume in the case of 
a sharp cutoff. Moreover, wide classes of smooth window 
functions give rise to the same asymptotic colored noise, 
which means that for choices of "good" window func- 
tions, physical quantities are independent of the window 
function. In fact, as was shown in [4J], the final cor- 
relations of the coarse-grained fields are independent of 
the window function if it is chosen to be 1) spherically 
symmetric and x-dependent only through the combina- 
tion |x|/_R (with R = [ea(i)7y(i)]~ , the length scale over 
which we want to coarse-grain), in such a way that Wu 
is constrained to have the form W^/f (x, i) ~ i?~^w(|x|/i?) 
when represented in configuration space and 2) such that 
the function w(s) is decreasing at least as ~ s~'° starting 



at s 



1. 



In what follows, however, we will be interested in iso- 
lating the effects coming from colored noises induced by 
self- and cross-correlations of the quantum fields them- 
selves, rather than colored noises induced by the choice 
of a smooth window function. Therefore, we leave Wn 
unspecified throughout this paper. 

With the question of how to define the window func- 
tion clarified, we can move on to describe how to remove 
the linearized equations satisfied by the quantum fields 
(12) (^^j T^j, on sub-Hubble scales from the equations of mo- 
tion for the full quantum fields $, ^. We can now see 
why a naive subtraction cannot give the right result. In- 
deed, because of the time-dependent comoving cutoff of 
the window function in Fourier space, as modes are evolv- 
ing in time, fc-modes keep on exiting the Hubble radius 
to join the mean, coarse-grained fields and in doing so 
escape the range of fc-modes where the quantum fields 
are defined. Therefore, upon subtracting the linearized 
equation obeyed by ^p from the full equation (|lip . we 
obtain: 



-Uip + m|(^ -f V;„rt,<i.((<5, x) = ^S'^-, 



y%t^{sp-.x)<^>i'-^ 



J^T rpcrt 

' 2 ,*** 



(V,X)0> 



2 ,**'!' 



(V', x)^l 



(13) 



Here, 5S^^ represents the effect of all the fc modes within 
8t which join the mean, coarse-grained field and give it 
an impulse (in the next section, we will see that this 
noise term arises from collectively treating all small scale 
modes inside a path integral; to give a rapidly-oscillating 
classical noise term). We can obtain an expression for 
(55^^ by substituting the expansion P^ in the second 
line of ([TT|) and noting that, since the (/i/c's are the lin- 
earized Fourier coefficients of $ for fc <C ali^ mode by 
mode they must satisfy the linearized Fourier transform 
of the total equation of motion (that is, 0fc -|- SHcpk — 



(fc2/a2-m| - K„t,$<i.)(?!'fc 



"Kort,**V^fe)- Wc find: 



'^^> " / (2.) 


3 




Hit))'"' 




^w ( '■ ^ 


X .. ( k \ - ' 

\ea{t)H{t)J 


-ikx 


" ^ lea(t)i?(i)y 






(14) 


^ke-^"'' = 





^flk + (j)\a_y. 


p-ikx 


(15) 



We can see how this expression relates to the qual- 
itative interpretation of SS^-^ given above as the im- 
pulse of every Fourier mode coming from the quantum 
field 0> when it joins the coarse-grained field by look- 
ing at the simplest example of window function, Wh = 
6 [k/{eaH) — 1]. Indeed, in this case, the time-derivative 
transforms the Wh it into a i5-function in the first two 
terms. This means that out of the integral over those two 
terms, only the mode crossing the window function at a 
given t contributes, and therefore these terms represent 
the kick to the mean field coming from the mode with 
k = ea{t)H{t). The last term includes the derivative of 
a ^-function, which might seem unsettling. However (as 
would be seen if we had been working with the Hamilto- 
nian first-order formalism rather than second-order equa- 
tions of motion Q), all these terms need to be included 
in order to represent the impulse to the coarse-grained 
field coming from both <j)k displacement and momentum. 

From the coarse-grained field point of view, we see that 
the quantum fluctuations escaping the Hubble radius act 
as a sustained noise source. Indeed, the statistical prop- 
erties of the (pk modes becoming larger than the coarse- 
graining region and constantly sourcing ip are that of a 
Gaussian-distributed (as the <f>kS satisfy the linearized 
Fourier space equation), zero-mean stochastic noise term 
^■^(x, t), with a variance which can be calculated from the 
variance of the quantum field. Note that, at this point, 
one should technically interpret this stochastic random 
noise as the piece of the quantum fiuctuations effectively 
acting as a classical field. As we will see in the next sec- 
tion, the piece remaining purely quantum will induce an 
extra term introducing a non-local mass renormalization 
of the coarse-grained field, as well as dissipation. How- 
ever, for the purpose of interest here where the slow-roll 
conditions arc satisfied, this term is negligible and we will 
not include it in the following. 

Defining 



?2 = / ll^^H [T^ijH 



(t) 
(t) 



Skflke 



-?k-x 



^k"-k'= 



it^ik-: 



Defining ^'''(x, i) in a similar fashion for the ^ field, we 
can obtain an identical expression for SS^^ . 

Note that here, it is strictly the time-dependence of the 
window function rather than mode-coupling (as would be 
the case in standard quantum Brownian motion) which 
gives rise to noise. Indeed, the noise only represents the 
continuous Hubblc-crossing of the modes, rather than 
some fundamental coupling between super- and sub- 
Hubble modes. This means that the properties of the 
noise terms generated in this way are such that they will 
not give rise to decoherence (in technical terms, a mass 
term in the equations of motions can only change the 
squeezing of a state, not induce and increase of its phase- 
space area or entropy, which means it cannot lead to its 
classicalization). See [13; ES] for an in-depth discussion 
of these points. 

In the case of the potential ([5]), we obtain the following 
equations of motion up to second order in the quantum 
fields: 



(16) 



it is straightforward to show that the noise term induced 
by the linear (/)> quantum modes exiting the Hubble ra- 
dius is given by: 



er - el 



-Kort:<I'**(</'>X)<^>V'> 



(18) 



°X + m|x + K.rt,*X 



-V 



***(V:X)^>'/'> 



ef 



?2 



iK 



p„t,***(¥':X)V'> 



5Kort,***('/':X)0> 



(19) 



SSd 



iHit + it-Q 



(17) 



The non-linear terms on the right-hand-side of both 
equations represent couplings between the coarse-grained 
fields and the small-scale quantum fluctuations. There- 
fore, one way of picturing the situation is to think about 
the coarse-grained fields as a system of long-ranged fluc- 
tuations immersed in a bath of small scale quantum fluc- 
tuations. The latter will back-react on the system and 
cause its decoherence, triggering the transition from long- 
range quantum mean fields to mean, coarse-grained clas- 
sical fields. In other words, the quantum mean fields 
become dressed by quantum fluctuations and they de- 
cohere through interactions with these quantum fluctua- 
tions around them, developing random classical fluctua- 
tions in the process [3J]. Note, however, that classical- 
ization is not only caused by inflation self-interactions. 
In fact, all non- linear couplings will contribute to this 
process (see, e.g., a treatment of classicalization from 
non-linear coup ling s to the gravitational tensor or vec- 
tor modes [3J, |43, or [i^l for an exploration of other 
possibilities). 

Equations (|18p and (|19p form a system of classical 
Langcvin equations, sourced by random gaussian noise 
terms (which are completely determined by their two- 
point correlation functions). They describe a stochastic 
process, which means that rather than solving for one re- 
alization of if and Xj we must solve for their probability 
distribution p(t, tp, x) over many realizations, through a 
Fokker-Planck equation [1J| . Expectation values of tunc- 



tionals of the stochastic fields, in particular their corre- 
lation functions, can then be calculated via: 



(F[(^(t,x),x(t,x)]) = / / duj^duJxp{t,uj^,uj^)F{uj^,uj^) , 

(20) 
where w^.^ are simply dummy variables for the fields. 

For completeness, the linearized equations of motion 
satisfied by the Fourier components of the quantum fields 
(j>k and ipk are: 



iH(l>k 



ipk + iHipk + 






Kcrt,** ) V'fc 



= -Vp„,,*<i,0fc 



(21) 



(22) 



From this, we can calculate the statistical properties 
of ^f-^. Solving the full system ^, ^ in the pres- 



ence of the classical coarse-grained fields (fT8| . (|19p is 
highly non-trivial, since not only is it a non-linear system 
with time-dependent effective mass for the two quantum 
fields, but the mean fields themselves are shifted by the 
presence of quantum fluctuations. In order to make the 
problem tractable, it is conventional to fix the quantum 
fields to be free, massless fields in exact de Sitter space, 
and to take their statistical properties to be s imp ly de- 
fined by the Bunch-Davis vacuum [H [H, [H [Ig. This 
approximation is justified at first order; however we are 
interested in studying effects such as the tilt of the spec- 
trum and the strength of back-reaction effects of small 
scale modes on large scale fluctuations; precisely the ef- 
fects that require a more precise knowledge of the details 
of the modes behavior coming from the quantum fields 
when they join the coarse-grained fields in order to be 
self-consistent. 

This is in contrast with most previous works who have 
been claiming to find non-perturbative solutions for cor- 
relations of the inflaton (e.g. |46l-l49f). by which they 
mean that they arc solving the classical coarse-grained 
system ^TE\i , ([T^ to arbitrary order in perturbation the- 
ory (in terms of the slow- roll parameters), in the presence 
of a noise ^ which is calculated from fixing the quantum 
fields to their exact Bunch-Davis values in de Sitter space. 
This approach of including higher orders in e in the classi- 
cal coarse-grained theory, whilst only using leading order 
approximations (in the slow-roll parameters) for the noise 
sourcing the classical perturbations is a priori inconsis- 
tent. In most cases of single field infiation, in particular 
for chaotic inflation, this nuance is just a technical point 
which doesn't really have measurable consequences, since 
the main source of non-gaussianities and the tilt in the 
spectrum are effects incorporated in the classical theory. 
Indeed, the tilt in {Sf, Sip) will then come from the fac- 
tor of H in ((T5)) and HH) , and the non-gaussianities from 
the derivative of the potential evaluated over the classical 
fields which is allowed to receive corrections proportional 



to the slow-roll parameters through the time-variation 
of the background ip as it rolls down the potential. As 
shown in, e.g., [4^, [SO], one then recovers agreement be- 
tween the "naive" stochastic method and the regular per- 
turbative approach. In particular, it has been shown [l6[ 
that for a scalar field with an arbitrary potential, one 
can in this way correctly capture all leading order loga- 
rithms [ln(a)]" of coincident correlators to all orders in 
perturbation theory. 

However, if two fields are present, situations might 
arise where it would be crucial that a correction to ^ 
is included. In particular, it would be the case in the sce- 
nario where one field is heavy or even develops a tachy- 
onic instability, as, for example, is the case in the final 
stages of hybrid infiation, where the value of the quantum 
fields correlators as they cross the window function de- 
viate significantly from their Bunch-Davis values. Even 
more importantly, the approach outlined above neglects 
the presence of back-reaction from the quantum fields 
onto the coarse-grained fields by setting the interactions 
between the two sectors of the theory to zero. In what 
follows, we will outline a tentative process to remedy to 
these flaws. 



D. Outline of the expansion strategy 

Our goal is to define a process which gives an increas- 
ingly precise result by including consistently every con- 
tribution at each order in the slow-roll parameters and 
in the fields. At zeroth order we start by solving the 
purely classical, homogeneous, noiseless equations of mo- 
tion under the slow-roll approximation. To zeroth or- 
der in slow-roll parameters, the small quantum fields are 
then approximated by their free, massless, de Sitter es- 
timates, as was done in the previous studies mentioned 
above. This makes finding the variance of the classical 
noise terms ^{'2 from the two point correlation function 
of (j)k and ipk possible, since in this case (PT|) and ([^ can 
be solved exactly in terms of the usual Hankel functions. 
Combined with the sharp cutoff choice for the window 
function, i.e. the theta function 9 [k/{eaH) — 1], we then 
find that the stochastic terms ^f 2 are exactly random, 
Gaussian, white noise terms, as expected. 

We can then substitute these values of the classical 
stochastic noise in ^TE\i and (|19p (keeping only terms of 
order lower or equal in slow- roll to that of the square root 
of the noise variance, the noise having zero mean by def- 
inition). Moreover, for the sake of consistency, since the 
stochastic noise incorporates exclusively the linearized 
quantum mode functions, we expand the coarse-grained 
equations to first order in the field perturbations in all 
other terms. This means all back-reaction terms from 
the small scale quantum fields onto the coarse-grained 
fields will be neglected at this order. The final result- 
ing coarse-grained evolution equation we find is then a 
Markovian process. We can solve this corrected equation 
for the statistical properties of the coarse-grained fields. 



In particular, we are interested in finding expressions for 
the mean value of Lp and x over many realizations, as well 
as their variance and more generally the whole probabil- 
ity density functions p{t,ip,x)- Effectively, we can think 
about this new classical background as being shifted from 
its zeroth order classical (or "bare") value by the small- 
scale quantum fluctuations constantly exiting the Hubble 
radius, freezing out, and joining the coarse-grained field. 
Qualitatively, it is easy to understand that, from the 
point of view of a quantum fluctuation deep inside the 
Hubble radius during inflation (more precisely, a fc-mode 
well above the cutoff scale) , the background it feels it is 
evolving in is not the "bare" , zeroth order background. 
Rather it is the collective background composed of the 
bare value plus all the quantum modes which already 
have frozen out and have become classical. This shifted 
background is exactly the one which the method out- 
lined above made precise and allowed us to solve for in 
the previous step of our method. Because of its stochas- 
tic nature, this solution is under the form of a probabil- 
ity distribution of the coarse-grained fields. This moti- 
vates the next step in our method, which is then to go 
back to (PT|) and ([2^ . and replace all occurrences of the 
coarse-grained fields by their average values, variances 
and higher momenta (this has been proposed, e.g. in 
[5l|,[52] and more recently in (SSj): 



f,X 






(23) 
(24) 



and solve the corrected linearized equations for the small 
quantum fluctuation fields, this time expanding to next- 
to-leading order in slow-roll (since we want to pick up 
the effects of the statistically shifted coarse-grained fields 
on the mass of the perturbations; in other words, we 
now solve for the full linearized mode functions). Qual- 
itatively, this corresponds to re-expanding to find the 
behaviour of quantum fluctuations around the adjusted 
background. From the results obtained for these cor- 
rected quantum fields, we can find the corresponding 
statistical properties of corrected stochastic noise term, 
which will now be colored rather than white. More pre- 
cisely, it will not be uniformly random, since it's maximal 
amplitude will be time-dependent, but it will again give 
rise to a Markovian process in the sense that the noise 
will still be memory less. 

From there, we can go back to the coarse-grained sys- 
tem ((T8| - (|19p . and reevaluate the probability distribution 
and the correlation functions of the classical fields in the 
presence of the corrected linearized noise, which now in- 
cludes non-zero correlations. To do this in a consistent 
manner, we now keep in the coarse-grained equations 

terms up to order \ {£,{2 ^it) i'^ slow-roll. This will 
allow us to find the spectrum of fluctuations (by solving 
for the variance of the classical coarse-grained associated 
field) predicted by stochastic infiation, including the tilt, 
in a consistent manner. 



This self-consistent implementation is expected to have 
important consequences, in cases where one or both of 
the fields are not test fields, in the sense that the noise 
amplitude depends on the stochastic value of the coarse- 
grained field itself - that is, when the amplitude of the 
noise is not (or cannot be approximated as being) time- 
independent^ . 

Two specific cases where we can expect stochastic ef- 
fects to play an interesting role when compared with re- 
sults previously obtained using the stochastic formalism 
are 1) in the neighborhood of an instability point in the 
potential, and 2) the case where the secondary field has a 
mass that is parametrically larger than the inflaton, such 
as in the case of quasi-single field inflation, or in hybrid 
infiation. In particular, in the case of hybrid infiation, we 
expect stochastic effects to play a preponderant role in 
driving the evolution of the background during the wa- 
terfall phase [la]. We derive the precise form of these 
effects in the second paper of this series [311 . 

Moreover, at this order in slow roll, we can now study 
the effect of back-reaction coming from second-order ef- 
fects in the quantum fields. In particular, we can es- 
timate the amplitude of the second-order terms in the 
quantum fields in (fT8|) - (|T9l) and see if back-reaction can 
have a significant effect on the coarse-grained fields. In 
the context of hybrid infiation, we are particularly inter- 
ested in studying this effect during the waterfall, when 
the ipy field becomes tachyonic and the V'fc modes enter 
a regime of tachyonic resonance. In this case, we expect 
the collective effect of the exponentially growing modes 
to leave an imprint in the background. If we find that 
the second-order terms become dominant over the lin- 
earized stochastic noise, this would be indicate a break- 
down of our perturbative approach in the quantum fields, 
which in turn can be used to rule out models. Alterna- 
tively, we can estimate the second-order contribution of 
the quantum field to the stochastic noise in the presence 
of the stochastically-shifted background, which can then 
be used as a criteria to measure the break down of per- 
turbation theory. We will come back to this study in a 
later paper. 



III. SCHWINGER-KELDYSH DERIVATION 

In order to define more precisely the idea of mean clas- 
sical fields arising from purely quantum fields, as well as 
to understand to what extend the primordial fiuctuations 



^In this case, one needs to assign a so called "prescription" 
a g [0, 1] to the Langevin equation in question, which sets at 
which point i + adt the noise must be calulatcd when the field is 
incremented between t and t+dt, when defining the Langevin dy- 
namics as a limit of a discrete stochastic pocess, and the Fokker- 
Planck equation associated with this Langevin equation will de- 
pend explicitly on a. We thank Vincent Vennin for pointing this 
out, and we plan to return to this point in a later study. 



are quantum and at which point they can be understood 
as providing a stochastic noise source to the background, 
we need to turn to the more technical formahsm of out- 
of-equihbrium quantum field theory. 



A. CPT formalism 

We once again start by splitting the full quantum fields 
$ and ^ into a long-wavelength and a short-wavelength 
pieces, with the cutoff between the two again defined at a 
scale around the Hubble radius. As before, we also define 
the small-wavelength part to be: 



'S'b,X,'/'>,V'>] = 5sys[i^,x] +'S'bath [(/>>, '0>] + 



S'int[</?,X,'/'>,'0; 



(26) 



0>(x,t) = /(f^Tyff(fc,t) 



ftkOke 



-ik-x 



'k^ke 



t ^zk-x 



V'>(x,i) = / ^WH{k,t) [tP^b^e-^^-^ + ri^bie 



tt „J'k-x 



(25) 



where Wh is still the same window function projecting 
out the long-wavelength modes which we will, again, pick 
to be a theta function in what follows. The mode func- 
tions 0k, ipk are chosen to be solutions of the linearized 
field equations, a point which we will make more specific 
later. 

The long- wavelength part of the fields will then simply 
be given by (/? = $ — </>> and x — ^ ~ '0>i again as 
before. Substituting this decomposition in the action for 
the scalar fields, we see that it splits into two distinct 
actions for each part of the fields, plus an interaction 
piece: 



J 



The fields ip and x represent the so-called system which 
we want to study, in the presence of 0> and ip^ which 
constitute the bath. In the following, our goal will be to 
perform an integration of the action over the bath degrees 
of freedom (the short wavelength fields), (/)> and ipy, in 
order to obtain an effective action Scfi for the system 
of long-wavelength modes, and finally take the classical 
limit. 

In order to obtain the in-in expectation values for out- 
of-equilibrium quantum field theory as opposed to the 
usual in-out transition amplitude, we proceed using the 
Closed Time Path (CPT) formalism [13, [H. In this 
framework, time integrals are made along closed contours 
starting at an initial time ti far in the past, evolving for- 
ward in time up to some time tg and then backward to 
the initial, distant time. The initial states are prepared in 
this distant past, point at which interactions and effects 
of curvature are assumed to vanish. This way, the states 
can be set to be their free vacuum state in Minkowski 
space, that is, the Bunch-Davis vacuum. We then split 
each field in the theory into two sets of fields, the first 
set for fields living on the forward piece of the contour, 
denoted by a ""'"", and the second one by "~" for fields 
living on the backward piece. This means $ is now de- 



composed into 4 scalar fields: (/?+, 
similarly for VP. 



•^ 



and 



and 



Since we want to integrate over the fields c/)>, 0>, -0>, and '0>, the object we are interested in is the reduced density 
matrix pr and its time-evolution operator jTi-. From the total density matrix p, we define the reduced density matrix 
Pr as: 



PriV^,X'^,V ,X ) 



d(j>+dtpi 



dcj)>dip-^p{Lp^,X^At^'^>'^ V 



,</>>, '0>)'5(0> 



6;)J(^+ 



i^>) (27) 



and is propagated in time from ti to to by the reduced evolution operator J^^'- 

/OO f-CG 

dvtdxt I dLp^ dxi Jriv'^ ,X'^ ^f] ,Xf ,tQ \ Lpj ,xt ,^i ,Xi ,ti)Pr{^t ,xt ,VT ^XT ,U) ■ 
-CX3 '^— OO 

(28) 
Here, the index / denotes the boundary values of the fields at the time ip- We would like to write J^r hi a path integral 
representation. This is in principle possible, however, in practice the explicit expression can be very complicated since 
Sr depends on the initial state. The crucial simplifying assumption we make is that at very early times ti, the system 
and the bath are not correlated, that is: 



p{t = ti) ^ Psys{ti) X /Obatli(ii) : 



(29) 



in which case the evolution operator Jr does not depend on the initial state of the system. This is a reasonable 
assumption in our case, since at sufficiently early times, the modes of the bath were deep inside the Hubble horizon 
and felt as if free in flat, Minkowski space. 

We can then write the reduced evolution operator in terms of a functional representation, and doing so we find an 
expression for the effective action for the system: 



± ,,± 



JriVf,Xf.tQ \y^^ .Xz ,U) 



Vtf^ 



4 



2?X cxp 



^± 



Vifi^ 



Vi 



4 



(^^ {S.ys[^+,X^] - ^sys[^-,X-]}) F[^^,x^] (30) 

(31) 



Xi 



^ PX*CXp<j^5eff[</?^,X^] 



where i^[v3 , X ] is known as the influence functional, produced from the action of the bath on the system. In general, 
it is a non-local, highly non-trivial object. Not only does it depend on the time history, but it also mixes the forward 
and backward histories along the time path in an irreducible manner. In our case of particular initial conditions, it is 
given by: 



F[^^,X^] 



+^d^p+^ / d0±_d^. 






exp ( ^{5'J^ath['/'>>'0>] - Sl^t\,[<P>A'>] +'5'intb+,X+,0>,^'>] -S'portb .X >'/'>>>]}) Pbath (<?!>>,, V'> J 



; exp 



-SiaW^,X^] 
n 



(32) 



Here we have used S'j^^th ^^ denote the action of the free bath degrees of freedom. Because of our choice of initial 
Bunch-Davis conditions, the integral of the initial density matrix over initial conditions of the bath fields can be 
assimilated to an overall constant (which we won't write in the following to avoid cluttering the notation). Note 
also that we have imposed the boundary conditions (j)y = (pZ, and ■(/'>, = V'>, at t = tg. In the last line, we 
have defined ^ia, which is called the influence action. Its relation to the effective action is easily found to be: 
^eff [(/?=*=, x^] = 5'sys[(/3+,x"'"] - SsysVp~ ,X~] + SiaW^^X"^] = 5'sys - •5'sys + S\A- We can write the infiuence functional 
as: 



exp<j -S'ia['^,x] 



dctyt^Mii 



(33) 



'/'>% 



'^ 2? -± / '^^ 2?V^±gi/rf''a;[(|0>A*0>-l'V^A,*0>) + (i0>Av.'A>+X^A^V^>)-"i4ath(0>,'A>)-"14„t(v 



where the time integral in the above spans from — oo to <0i the turn-around time in the CTP contour^, and where we 
have introduced the matrix notation for simplicity: 



V-; 






V 



A. 



A0 
-A. 



A„ 



Av, 
-Av, 



X+ 

x> 



(34) 



with the superscript T meaning transpose, and we have also introduced A^ and A^, which are defined to be the 
integration kernels for the free scalar fields: 



A^ = ~a^{t) 



d^ + 3Hdt 



Hi) 



A^, 



-a'{t) 



d1 + 3Hdt 



a\t) 



(35) 



^However, since we impose the boundary conditions ^^(x, io) = 
$~(x,to) and \l/+(x,io) = ^~(x,4o)) we could easily extend the 
upper bound of integration to future infinity, as the forward and 
backward parts of the contour cancel exactly beyond to- Indeed, 
to being the time at with we are evaluating a particular operator, 
causality of the theory requires this cancellation. 
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In the first two terms in the exponential we put all the terms which are bilinear in the fields, i.e. the terms over 
which we can compute the path integral exactly to obtain the leading order influence action, and in the last two terms 
we put all non-linear couplings, for which in general we need to do a perturbative expansion in the coupling constants 
in order to perform the path integral over the bath degrees of freedom. Note that V with any subscript are the vector 
form of terms which were part of V^pcrt as was defined in the previous section, as all potential terms in V involve more 
than two-fields interactions. In the case of the generic potential (O, this will generate both loop corrections coming 
from the self-interactions of the baths degrees of freedom, represented above by the term Vbath, and loop corrections 
due to the interactions of the bath with the system, represented by the term V^nt- 

B. Schwinger-Keldysh derivation of the Langevin equations 

We want to start by calculating the leading order influence action 5*1^ in order to obtain the propagators (or 
equations of motion) for the system degrees of freedom once the leading interactions between bath and system have 
been integrated out (note that at zeroth order, the influence functional trivially only adds a multiplicative factor to 
the system path integral). Therefore, at this order, we neglect the interacting potential both in the bath and at the 
bath-system level beyond the terms linear in the interactions, and only perform the exact path integral over bilinear 
terms in p3p . as done in detail in Appendix A. That is, only taking into account linear and quadratic terms in the 
influence action, we obtain a gaussian integral which can be performed to yield: 



S[^j^ = ^ f d^xd^x'ipq{x)Re[n^{x,x')]ipq{x') - | f d*xd^x'9{t 



t')ipqlm [U^{x, x')] fc + ix^ V') 



(36) 



where: 



n0(.T,.T') = 



(27r)3 



a\t)[PtMt)]e-'''-''a-'{t')[Pt,ri,{t')] 



J./M „ik-x' 



(37) 



with Pt defined as: 



Pf 



Wnit) + 3HWHit) + 2WHit)dt 



(38) 



and where we have moved to the so-called Keldysh basis and have defined the classical and quantum fields: ipc, Xc 
and ipq, Xq, respectively, by [l3|: 



fq 



<p +'P 



^ 



+ 



^~ 



Xc 
Xq 



(39) 



r 



In this manner, we reproduce the result of [56|. Note 
however that here the mode functions appearing in 
T\^^^{x,x') are technically calculated from solutions to 
the full massive mode equations of motion, since we con- 
sidered the masses to be part of the tree-level propa- 
gators as opposed to perturbation. This will be of im- 
portance later, when we include interaction and obtain 
system fields as part of the mode function equation for 
the bath degrees of freedom, since in that case the mode 
function appearing in similar expressions will again be 
calculated from the full mode function equation of mo- 
tion, with the actual system fields (as a justification for 
this, recall that all the operations done here are inside a 
path integral for the system fields), rather than on some 
fixed artificial background. 

The leading order of the influence action which we 
defined in ^ by F[if^,x^] = exp ^s'l^ therefore 
splits into a real and imaginary part. As pointed out 



for example in |57l . ISSj , they respectively represent noise 
and dissipation in the system fields. Indeed, the ker- 
nels Im[n0(x,x')] and Im[n.^(a;,a;')] are known as the 
dissipation kernels for each of the classical fields. The 
motivation for the name comes from the fact that their 
non-symmetric part each introduces an extra term in the 
equation of motion of each classical field when taking the 
variational derivative of the effective action, SScs/Sipq 
and 6Soff/5xq- These terms are non-local and propor- 
tional to ipc and Xci respectively, and can therefore be 
assimilated into friction, or in other words dissipation. 
However, as long as we are interested in studying the 
equations of motion while the slow-roll approximation is 
valid, these friction terms will be negligible in compari- 
son to the Hubble friction terms, related to the Hubble 
expansion (as shown, for example, in (56|, |59|). 

Similarly, the symmetric part of the kernels 
Im[H0(a;, a;')] and Im[H^(a;, x')] each give rise to a mass- 
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rcnormalisation term which does not contribute to the 
mixing of the forward and backward histories 1601 . These 
terms can also be shown to be neghgible [5a, [S^l ■ There- 
fore, in the foUowing, we will neglect the contributions 
to the effective action Scs proportional to Im[n0(x,x')] 
and Im[n^(a;,a;')] altogether^. 

On the other hand, the terms containing the real part 
of n^ and n^, in p6p give an imaginary part to the ef- 
fective action. As such, this term cannot be properly 
interpreted as an ordinary part of the action. That is, 
as detailed in Appendix A, this terms can be reinter- 
preted as a result of a weighted average over the con- 
figurations of stochastical noise terms representing the 
sub-Hubble quantum fluctuations which couple to (p and 
X- To make this explicit, it is possible to re-introduce 
these as classical noise terms with the appropriately cho- 
sen probability distribution [29|, as is customary to do 
in out-of-equilibrium condensed matter physics or when 
studying Brownian motion jl7| . We will do just this in 
the following. 

Before we do so, we first comment on how this in- 
terpretation of the real and imaginary part of the in- 
fiucncc action can be linked with known results from 
out-of-equilibrium statistical mechanics. As in the the- 
ory of quantum Brownian motion, the presence of quan- 



tum fluctuations gives rise to stochastic effects in the 
system through two distinct effects: the first one is 
a rapidly-oscillating and stochastically distributed noise 
term in the systems fields equations of motion; the sec- 
ond one proportional to the velocity of the system fields 
is a friction term describing how fast the system ap- 
proaches equilibrium from an out-of-equilibrium config- 
uration. The latter depends on the correlation of the 
stochastic force, and is thus not independent from the 
noise term. In out-of-equilibrium statistical mechanics, 
this relation is known as the fluctuation-dissipation theo- 
rem. It is evident here because the dissipation and noise 
terms are related through the real and imaginary parts 
of the same kernel. Typically, in non-equilibrium stati- 
cal mechanics, this theorem further states that the ratio 
of the transport coefficient and the microscopic random 
fluctuations is proportional to the temperature of the sys- 
tem. One could therefore ask whether in our case the 
ratio of the two terms we flnd here can be related to the 
Hawking temperature H/2Tr of de Sitter space. As shown 
in [5g, it is indeed possible to define a 'temperature' of 
order of the Hawking temperature, but (in the case of 
quasi-de Sitter space) which depends on v and e, and 
therefore is a property specific to the matter-spacetime 
interaction. 



J 



Here, in order to be able to properly interpret the imaginary part of the infiuence functional as noise, we introduce 
two real classical random fields per field in the system ^{ , ^2 ^-nd £,{' , ^21 each obeying the Gaussian probability 
distribution: 



?1 ''52 



exp^-i / d^xd^x'[^t'^{x),^p^'ix)]A-\x,x') 



^P^'ix') 
l^P^ix') 



where we have defined, letting r = |x — x'| and x = (t, x). 



Hjr) 
m{T) 



H{r') 



Re 



where 



M^\kT,kT') = 



Mmut') 



q{^kT)Mt)Kit') 



Mt)Kit')q (-fcr') q (-fcr) MtWkit^q {-kr') 



(40) 



yVi\kTM') , (41) 



(42) 



and similarly for the matrices AV(x,a;') and MV(A:t, fcr'). The functions q4,^^{—kT) represent the fc-dependent (and 
time-dependent only through the combination k/aH) part of the logarithm derivative of the mode function of the 
quantum fields, as defined in (|A lOp . Also, note that, again, the mode functions appearing in MV(fcT, fcr') are 
calculated from the full mode functions equations of motion. 

After some manipulations (see Appendix A), we can now use the infiuence functional to rewrite the leading-order 
effective action. To do so, we recall that we are interested in taking the classical limit of the effective action. As 
mentioned above, the fields tpq and Xq embody the quantum nature of the difference between the configurations on 
the forward and backward pieces of the time contour. In order to explore the classical limit of the effective action, we 
proceed to the field rescaling (fq, Xg ^ ^Vqi ^Xg [III- With this redefinition, a first order expansion of the effective 
action in hipq, hxq corresponds to the /i- independent, or classical, sector of the action: 



^In particular, |5a | argued these terms to be 0(e'^), whereas the 
Re[n0 ,^(3;, x')] terms we keep here are at most O(e^). 
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Vip^ 



4 



Vx cxp 



* c.(l) 



+x, [(□ - ml) xc - Ko.,* (^c, xc)] + ^, [p0(t)ef + ^2^] - ^,ef + X, [p^mt + el] - x,ef }) 



, (43) 



r 



where Pif,^^ are defined as the /s-independent part of the 
derivative of the mode function of the quantum fields. 
Note that the separation of tlie infiuence action ([55)1 we 
have just done between a weiglit part ^'[■^i,^] and a 
pure action part contributing to the effective action S'^g 

is unique [23 ■ The fields SX ^"^^ it can then easily 
be interpreted as classical random fields, since they obey 
Gaussian statistics and their weight function has a real 
exponent. 

The equations of motion for the classical fields ipc and 
Xc in the classical limit are then given by: 



,55'' 



(1) 



SLPq 



ss. 



(1) 
off 



^,=0 



SXq 



(44) 



x,=o 



Note that we evaluate these for iy9g = Xg = since our 
choice of boundary conditions at i = ig enforce that the 
actual time axis is unique. Relabeling to ^ t for simplic- 
ity, we therefore obtain: 



(-n + m|)(^c + V-s-i^c, Xc) 



t<t> li'P 



p4m +el +ef + 3i/ef,(45) 



(-n + mDxc + V^,{(pc, Xc) 



pAmr HI' +er+3ffer.(46) 

Here, the Gaussian noises are defined to have zero mean, 
and therefore their statistical properties are completely 
determined by the matrices A^-* {x,x') and A^f {x,x'). In 
particular there can (and will) be cross-correlations bc- 



t'P'i' 



tween y '^ and ^f , and their self- and cross-correlation 
functions will depend both on the choice of window func- 
tion and the level of precision we use to specify the mode 
functions (and their derivative through the g^^^ i~f^T) 
functions introduced in (jH])) at Hubble-crossing. This 
is of crucial importance for the implementation of the 
solution method we outlined at the end of the section HIl 
Comparing the result we just obtained in dJS]) and pS)) 
with p^ and ([T^. we find that the stochastic terms are 
exactly consistent, since if we go back to the way we 
defined Ci ou ^^'^ it old i^ section HTJ we have: 



it 



old 



d^k . 



'2 old 



(27r)3 
(2^ 



Wncl^ke 



:WH(t>\^e 



-■ikx 



-■ikx 



(47) 
(48) 



and similarly for ^2 ■ However, recall that, along the way 
of the derivation of the stochastic equations of motion we 
presented in Appendix A, we assumed the time derivative 
of the mode functions of the quantum fields, (pk and "fAfc, 
could be written as: 

•^k = [p,p{t) + q^ (-fcr)] 0k. (49) 

We therefore find that Cf old ^^"^ ^2 old "^^^^ ^^ spVii into: 



^1 old '>! new SI 1 

e2^id = -p.mt-it 



(50) 
(51) 



and similarly for Cioid 



and (,. 



2 old' 



second term of ^2 -i-^ ^^ ^^^^ form, 



^2 old 



where, to obtain the 
we have absorbed 
the g^^^ i~kT) functions in the matrices AV(a;,a;') and 
AV(a;,a;'), as we have done above, and where we have 
used that Wr ~ Hkr (l + H/H'A dkrWn- 

We therefore obtain the perfect consistency between 
the two methods, up to linear order in the bath fields. 
In addition, it is very easy to extend the above method 
using the usual tools from field theory to include higher 
order (in the bath degrees of freedom) corrections to ps)) 
and (|46p . and understand effects such as back-reaction. 

We now want to understand how this calculation can 
be done in the presence of a time-dependent window func- 
tion (see e.g. [23, |30, |6l| for a calculation of loop correc- 
tions up to second order in the context of a single scalar 
field with a power-law potential in Minkowski space and 
a time-independent window function) , and see how these 
corrections give rise to the back-reaction terms we found 
in (dll) and ^. 



IV. ONE-LOOP CORRECTIONS 

A. Setup 

In order to illustrate how loop corrections can be un- 
derstood and calculated, it will be useful to assume a 
particular form of interaction potential between bath and 
system degrees of freedom, Vint(<^,X, <^>, V'>)- In partic- 
ular, this boils down to choosing an interacting potential 
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for the fields $ and ^, which we choose a renormaliz- 
able, Z2-syninietric interacting potential in the symmetry 
phase, as follows: 



V^pert(<I>,*) = -* 



-$^\I/^ 



(52) 



It will also be interesting to consider the case of a mass 
hierarchy between the two fields: to$ <?C mq,. 

The first thing to do here is then to split the full fields 
into bath and system degrees of freedom, as we have done 
above, and to expand perturbatively in the bath fields, 
with the goal of performing their path integral order by 
order in a loop expansion, in a fashion very similar to 
what is done in Wilsonian renormalization [621 163| . 

But, before we do so, at this point we need to define 
the free action for the bath degrees of freedom. Indeed, 
we are facing something which looks like an ambiguity 
here, because once ([52)) is expanded in system and bath 
fields, (f, X s-nd (/)>, ■(/'>j respectively, it gives rise to in- 
teraction terms which are bilinear in the bath degrees of 
freedom, but contains further powers of the system fields. 



such as e.g. a jX^i^> term. We need to decide whether 
these terms shall be included in the propagators of the 
free bath fields as time-dependent masses, or treated as 
perturbations, as explained in [30| . Happily, it is possible 
to show that these two methods are equivalent, as will be 
explained later and in Appendix C. For now, let us sim- 
ply argue that, physically, quantum fluctuations which 
are part of the bath (and so sub-Hubble) , are evolving 
on a background which they feel to be the zero mode 
plus all modes which have already crossed the Hubble 
radius and which we are now collectively treating as the 
coarse grained fields, even at tree level. It is therefore 
intuitive to include terms which are bilinear in the bath 
fields such as jx'^i^> &s time-dependent masses in the 
free propagator of the bath degrees of freedom. 

Finally, let us also stress that one should really use 
the actual system field in the bath propagator, since we 
are doing the path integral over the bath fields inside 
the path integral over the system fields. Another way to 
see that this is indeed the case would be to recast this 
problem in the language of the 2PI formalism [gj, Ifia] , 
which we leave to a further study. 



J 



With these issues clarified, and going back for a moment to the +/— notation for the fields rather than the Keldysh 
basis for the classical/quantum notation, we can now move on to writing down the free bath and interacting action 
corresponding to ([5^ for the -I- and — parts of the contour: 



cO+ _ qO- I c± 
'-'bath '-'bath ^ '-'int 



^Ht) 



2g'^+x^<^ti't+9' {^^y<Pt m'+9' ix^Yi't m' + ^x^ii'tr 



{+ 



o[(^±)^(0±^^ 



'•^^^ '-^'4r] 



-} 

(53) 



where we have introduced a matrix notation similar to the one introduced in (1341). with: 



a(1),(2) _ 



A 



+ (1).(2) 





-A 





-(1),(2) 



<P,i> 



(54) 



and the integration kernels to first and second order in the bath fields are defined to be: 



^±.(l)^^±,(2)^^3(i)rn_^|_^^2(^±)2 



A 



±.(1) 



a3(i)[n-m|-g2(^±)2_ A(^±)3- 



A^ 



±,(2) 



a'it) a^ml-gH^^y-^ix^y 



(55) 



r 



B. Unperturbed influence functional and free 
propagators 



In order to obtain the loop expansion and calculate 
the effective potential for the system degrees of freedom 
once the bath fields have been integrated out, we intro- 
duce four currents, one for each field on the forward part 



of the time contour, and one per field on the backward 
part, as in the standard procedure [20, |61|. We then pro- 
ceed as outlined in Appendix B to define a perturbative 
expansion of the corresponding influence functional. 



Evaluating the first term in this perturbative series 
yields the unperturbed influence functional F^^'[J^, J^], 
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as defined in (jB 2p . which can be explicitly calculated: 

-1/2 



These propagators are represented below: 



F^^\j4>Ji,] 



exp 



InDet 



{^ 



(56) 



-^ I d^xd\'.rAU'^ 



-1 _ 

J0 + (0 ^ IP) 



Note that here, ( A j, I can be defined in terms of 

Green's functions on the contour, in analogy to (jA 3p . 
Also, because of our definition of the free propagators 
for the bath degrees of freedom, the determinant terms 
in the above arc not simply an extra multiplicative con- 
stant added to the influence functional, and cannot be 
neglected as is sometimes assumed. 

From here, it is easy to verify that the one point func- 
tions of the bath fields gives zero, that is, 

(0±(x))o = (V^(.t))o = O. (57) 

Furthermore, we can compute the two-point functions: 

(0>(a;)0?:(2/))o (58) 

T(<^+(x)(^+(y)> (0;(y)(^+(x)) 

(0;(x)(^+(y)) T(0;(x)0;(y)) 



and similarly for ■(/;> . As expected, we also find 
(,^>(a;)^>(2/))=0. 

Here to calculate these correlators we need to ex- 
pand the fields in Fourier space in terms of creation and 
annihilation operators. The mode functions appearing 
then need to be solutions to the equations of motion 

Aj''(fc)0k - 0, Aj)(fc)^k = 0, where A^Ufc) are de- 
fined as the Fourier transform of the kernels defined in 

dMli-dSS]). 

Now, since we are interested in writing down a di- 
agrammatic expansion, the above correlators are inter- 
preted as the bath propagators iftA0(a;, y). More specifi- 
cally, for fields time-ordered along the forward CPT con- 
tour, they are equivalent to the standard time-ordered 
Feynman propagator, ihAp{x,y), while for anti-time- 
ordcrcd fields along the backward part of the contour 
they correspond to the anti-time ordered Dyson prop- 
agator ihA]j{x,y). The upper off-diagonal element of 
((58)) represents bath fields propagating from the -I- to 
the — pieces of the contour, and therefore has x'^ al- 
ways 'earlier' than y^. As such, it is proportional to the 
absolutely-ordered negative-frequency Wightman propa- 
gator ihA^(x,y). Similarly, the lower off-diagonal ele- 
ment of ([58]) represents bath fields propagating from the 
— to the -|- pieces of the contour, and analogously is pro- 
portional to the positive-frequency Wightman propaga- 
tor ihA^{x,y). Therefore, we can identify: 



(0>(a;)0>(y)) = ihA^{x,y) = 

^<f,F{x,y) -A^^<{x,y) 
-A^^y{x,y) A^,D{x,y) 



(59) 



+ 



+ + 



:« <*: 



+ + 



; O CP ; 



We use double lines to represent bath degree of free- 
dom, and reserve the use of single lines for system fields. 
Dashed lines (either double or single) are used to repre- 
sent $ fields, while plain lines are used to represent ^ 
fields. Moreover, in the following we will use solid dots 
to represent vertices on the forward piece of the contour, 
while open dots will be used to represent vertices on the 
backward piece of the contour. Similarly, the signs '-I-' 
and '— ' indicate if the field represented sits on the for- 
ward or backward piece of the contour. To simplify the 
notation, for bath fields we will only indicate '-|-' or '— ' 
on the appropriate line, while for system external legs we 
write explicitly which field the line belong to. 

From here, converting the bath fields propagators to 
the Keldysh basis, where the fields are expressed in terms 
of classical and quantum fields, we find: 



ihA^{x,y) = 

4Re \^^{x, y)\ 2ie{y'' - xO)Im [n^{y, x)] 

2iQ{x°-y'')Im[n4x,y)] 



where 



^4'i^,y) = 






(60) 



Tyi,<0k(x°)(^^(y°)e"^'^-(''-'') (61) 



and similarly for ■!/)>, with 0.0 (x,?/) analogously defined. 



C. One- loop calculation 

Having found an expression for the bath fields propa- 
gators, we can move on to computing the remaining one- 
loop contributions. Note that diagrams involving vertices 
with only bath degrees of freedom will give their first 
contribution at two loops, and so we conclude that terms 
involving Vbath(0>, ''/'>) "^iH give no contribution at one 
loop. The same goes for terms C (■(/'>, V'> </'>)■ From ([55)) 
we see that the first contribution to the loop expansion 
(see ()B 6[) for the explicit details of the expansion) will 
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come from terms involving two powers of the interacting 
action: 



F[ip^,X^ 



""n 



^(^.M.-,x-,0i,^>-);„ 



-2^(">^nt(¥' >x ,4>>,i'>) 



-^ (Mnt(¥'^,X^,0>,'0>)Vint(<^ ,X >>,'0>) 



coming from the fohowing terms in the action: 



(62) 
(63) 



We now work out the impact of these two terms one af- 
ter the other. As before, one would argue that the two 
terms in (jHH) , which are linear terms in the bath degrees 
of freedom, vanish, since the elements of the pairs {if, (j)} 
and {Xj "0} s-rc defined to be orthogonal to each other. 
However, in our case, the kinetic term in A^ < acts on 



the window function in (^> and ?/>>, making the contri- 
bution from these terms no n- vanishing. This is the exact 
effect we calculated throughout section lllll and with no 
surprise we find their contribution to be 



2h 



d^xd^ 



a^{t)(plDxa3A^{x,x')a3Dx'(fx'a^it') 
+{^ -^X, 4>^^}], (64) 



where ua is the third Pauli matrix. Transforming to the 
Kcldysh basis, and proceeding through a series of manip- 
ulations similar to the ones leading to (jA 12[) in Appendix 
A, we recover (|36p . Therefore, reorganizing the perturba- 
tive expansion, we find that (j62|) gives rise to a correction 
to the effective action identical to the one calculated in 
(|33|), with ^{'2 interpreted as random classical Gaussian 
noises with variance given by the probability distribution 
(pn|) . The only difference here is that the mode functions 
used to calculate the MV (fcr, kr') matrices given by (^^ 
must be solutions to the equations of motion including 
interactions, K^^\k)(j)i, = 0, A^^^(fc)V'k = 0. 



Then, equation (j63p leads to a contribution at one loop given by the diagrams: 



^\. 



\,.'^"<t./ 



+ /^^ <^t + _ 



• f V, 



\,y~<t./ 



- _ + / 



if' if^ 



ff 



\/''~<'../ 



yvK/ .y^.-\-^=A..\-AA 



'X X 



X X 



'X X 



Evaluating these graphs, we obtain 

■,4 



^^, Jdxdxa {t) [^ {x)x {X), ^ {X)X (x)J (^ _(^-^+)(^-^+) f{<t>Utm^ti^i) 



^+{x')x+{x') 
ip {x')x {x') 



3r^'\ 



a^it 



Rotating to the Keldysh basis, we obtain, up to second order in the quantum fields (terms of third and fourth order 
in Xq and ipq will also receive contributions from two and higher loops, and give higher contributions to the one-loop 
Langevin equation, we therefore neglect them in the following), the following contribution to the influence action: 



(65) 



ggl-loop 



(66) 



^^xd^-'"^' 

-,4 



2g^- I d^xd'^x'a^t) {^,Xx)xq{,x) + ifq{x)xc{.x)) Re [%{x, x')n^{x, x')] {<f,{x')xq{.x') + ifq{x')xc{x')) a^t') 



d'^xd^x'a%t)ip^{x)xc{x)e{t - t')lm [f]^(a;, x')n^{x, x')] {ipc{x')xq{x') + ipg{x')xc{x')) a^{t 



As before, we obtain both a real and an imaginary contribution to the influence action, which we interpret as a 
a noise and dissipation term, respectively. We again see a manifestation of the fluctuation-dissipation theorem. As 
before, the dissipation term will be negligible in the context of slow-roll inflation where most of the dissipative effects 
come from the Hubble friction term in the equations of motion. Also, again as before, in order to properly exhibit 
the noise as so, we introduce the classical random noise terms ^3 and ^3 defined by the probability distribution: 



V 



?3'S3 



exp 



d^xd^x'{5X{x),it(x)\Q-\x,x') 



C|(^') 
Vit{x') 



(67) 
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where, if we assume that the external momenta of the coarse-grained fields are small compared to the internal momenta 
running in the loop, the B matrix can be simplified as follow: 



/al\ 77 
—^e [Wh (kr) Wh {kT')f Re [N*'^'(fcT, fcr')] (^^(x - x') + c.t. , 
ZTT 



(68) 



with 



W^^{kT,kT') = 0kW'/'k(iV-k(i)^-k(i') ( \ \ 



(69) 



r 



Although IR finite thanks to the presence of the coarse- 
graining window functions, the integrals in (j68p will be 
UV divergent once we plug in an explicit mode function 
with Bunch-Davis initial conditions (as opposed to the 
integral in (HJ), which was both UV and IR finite due to 
the presence of derivative of window functions, effectively 
mimicking the effect of delta functions). However these 
UV divergences are easily dealt with using the standard 
procedure by introducing a physical UV cutoff A and a 
local counter term (which we denote c.t. in the above), as 
in [6g| . This translates into a UV cutoff in the comoving 
momentum k at aA, where the scale factor a is evaluated 
at the smallest of the two times t, t' . The latter must then 
also be regulated with the same physical cutoff, which is 
most easily done by rewriting the time integrals as: 



dtdt' 



dt / dt' 



dt 



dt' (70) 



The local counter term will then cancel any power law 
divergence, while replacing logarithmic divergences by 
terms proportional to In I — j + C, where /x is the renor- 
malization scale, and C is a constant. 



The contribution to the influence action coming 
from the term containing the real part of the kernel 
5l0(a;, a;')ri^(x, x') in (j66p then becomes: 



^gl-loop 



(71) 



2g^ / d'xa\t) it{x), S.t{x) 



Vc{x)Xq{x) 
Xc{x)Lpq{x) 



D. Contributions of the InDet terms at one loop 



To write down the one-loop Langevin equations, we 
only have left to deal with the two In Det(A\_ ^)~^^^ terms 
in (j56p . as they contain resumed powers of the system 
fields which need to be differentiated in order to obtain 
the right equations of motion. Calculating its influence to 
first order as done in Appendix C, we obtain a non-local 
mass renormalization term, a new noise term, as well as a 
new dissipation term. Once rotated to the Keldysh basis, 
we obtain for the (p and x fields, keeping only corrections 
up to one-loop order: 



\\nJ^ei{Kf) 



= -^ \ -y / d^xa^{t)xcXqF<p<p{x,x) 



h 

23^ 



+g^^ / d'^xd'^x'a^{t)xc{x)xq{x)Yic[F^^,^^{x,x')]xc{x')xq{x')a^{t') 



d^xd^x'a^{t)xc{x)Xc{x)e{t - t')lm [F^^,^^{x, x')] Xc{x')Xq{x')a^t') \ , 



(72) 



ilnDet(A(,^)) 



A q^ 

d^xa-^{t) [ -XcXq + l^^cVq ) F^,^{x,x) 



+ ^ I d^xd^x'a^t) 



g'^ipc{x)lfq{x), -Xc{x)Xq{x) 



d'^xd'^x'a^{t)9{t - t') 



A 



g ipcix)ipcix), -Xcix)Xcix) 



Re[F^^_^^(x,x')] 
Im[F^,^_^,^(a;,a;')] 



g'^ipcix')(pg{x') 



\Xc{x')xq{x') 



g^ipc{x')ipq{x') 

^Xc{x')Xq{x') 



(73) 



a3(t' 
a^t') 



where the different _F's above are related to different powers of the propagators of the free, massive, bath theory. In 
our case, for the purpose of giving an example, we can write them using the massless mode function for the inflaton 
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field since we assume it to be light, and the fuU, massive mode function for the spectator field. Their regulated 
expression can then be explicitly written as: 



FMx,x) = f ^P(WhY^{-rr\H^{-kr)f + c.t., 



dk ,„^f n2 -ff , ^ 



c.t. 



(74) 
(75) 



Ftj)tj),tj)tj){x^x') = 



ah 



^e [Wh (kr) Wh {kT')f Re [0^{kT, kr')] S^^ - x') + c.t. , 



■ ipip,ipip 



{x,x') 



"^ dk 



^k^ [Wh (kr) Wh {kr')]' Re [o;;^ (fcr, fcr')] ^^(x - x') + c.t. 



where 'c.i.' denote the contribution of a counter term, and where: 



(76) 
(77) 



OJkT, kr') = ^ ri + fcVr' + ik{T - r')] ' e'^^'^t"-"') 
El^irrr[HA-kT)H:{-kr')f (\ J 



O^^" (fcr, fcr') 



(78) 
(79) 



Here, we use f to denote the usual mass-correction factor in the Hankel functions: i^ = \/9/4 — '^i^/i/^. 
Introducing the extra classical gaussian noise terms ^4 , ^4 and ^5, defined by the probability distributions: 



e4> c4' 
'54 ''54 



exp^-- / d^xd^x'[^t{x),^tix)]F^l^^{x,x') 






V [^5] = cxp 



-i J d^xd^x'U^)F-l^^{x, x')i,{x') 



(80) 



and neglecting the dissipation terms in ([75)) and (17^ . we can then rewrite the two lnDet(A^ .,.) ^/^ terms contribution 
to the influence action up to one loop as: 



rf,,^) 



ss 



InDct 
lA 



d xa (t) <- { -XcXq +9''Vc^q F^^{x,x) - g XcXqF4,4,{x,x) +2 



(:<P fV 
S4 ^ S4 



g'^ifc{x')^q{x') 

^Xc{x')Xq{x') 



(81) 



+ '2'9^£.5Xc{x)Xq{x) 



Therefore, bringing all the pieces together, we can collect the above contribution to the influence action together with 
what we calculated from the one-loop interaction in (j7ip and add all these interacting one-loop corrections to the 
free, leading order result we have calculated in (|43p . Since we are, as before, interested in exploring the classical limit 
of the effective action, we once again proceed to the field rescaling ^Pg, Xg — > f^Lpq, hxq. We again want to perform a 
first order expansion of the influence action in hipq, hxq, which corresponds to the ft-independent, or classical, sector 
of the action: 



Vip"^ 



Vx exp 



«o(l) 



0/0 



= Pv'^'^Px''^ / 'Deii.A'^i^'p ii.i. 



P,i> i:4>,''P 



S3, 4' S3, 4 



^Ks 



(82) 



X exp f i / d'^xa^{t) [ipq [{U - m%) ipc - "Import,* (V'c Xc)] + Xq [(° - ml) Xc - "V^pcrt,* (<<5c, Xc)] 
'.mt + it] - ^qit + Xq Umt + it] Xqit + 2.9^ ^ (^) > ^^ (^) 



+ '^q 



-^XcXq + 9 (Pc'^q) F.^^{x,x)- g XcXqF<l>4>ix,x)+'2 



S4' ?4 



'^^c{x')Lpq{x') 
\Xc{x')Xq{x') 



ipc{x)Xq{x) 
Xc{x)ipq{x) 



+ '^9'^£.5Xc{x)Xq{x) 
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The equations of motion for the classical fields ipc and Xc are then obtained via the saddle-point conditions (|44)) : 



(-D + m|)^, + V;$(^„ Xc) = P^mt + ^t +it + ^H^f + 2g^£,t{x)^,{x) + 2g\t^,{x) 

'^'X.c-fipijiyx^ X) , 

{-D + ml)xc + V^{^c,Xc) = pAt)^t+^t+it+3mt + 29'^ti^)Xc{x) 

-g^(PcF^4,{x,x) - g^XcF4,4,{x,x) + X^txc{x) + 2g^^5Xc{x), 



(83) 



(84) 



where the mode functions used to find the variance of the various noise terms in the above are the one solving the 
equations of motion: 



O-ml-g^x^)" 



n-ml-g^^^f-^{x' 



Mt) = 0. 



(85) 



r 



This is the main result of the paper. We clearly see that 
the extra terms derived at one loop are actually introduc- 
ing back-reaction effects in the sense that they couple two 
bath fields with a (classical) system degree of freedom as 
a source in the classical fields effective equation of mo- 
tion. When interested in the problem of back-reaction, 
the terms appearing in the above Langevin equations will 
be the leading ones in powers of the coupling constants 
of the theory embodying those mode coupling effects. 

Linking back to our heuristic derivation from section II, 
these are the leading contributions to the extra terms in 
the Langevin equations we found in (fT8|) and ([T9l) . More 
specifically, these are the precise terms which were second 
order in the quantum fields and which we wrote explicitly. 
This means that these second-order mode coupling terms 
in the quantum fields can in fact be rewritten as a sum of 
three different noise terms for the x fi^ld and two different 
noise terms for the ip field (closely linked to the fact that 
there are three really distinct diagrams at second order 
in the couplings contributing to the x effective potential, 
while there are two for the ip effective potential). 

It is interesting to see that in the heuristic case, 
the method used could not capture the contributions 
to the equations of motions coming from the mass- 
renormalization terms in the effective potential of the sys- 
tem fields (which are the terms on the RHS of the above 
without noise terms), as well as the dispersion terms, 
which ended up being negligible in the slow-roll approx- 
imation and which we neglected in the above. Note that 
this was also the case when we evaluated the free in- 
fluence functional in section III, where we obtain a real 
term for each system field giving rise to dispersion but 
decided to neglect it because their contribution was again 
negligible under the slow-roll assumption. The heuristic 
approach did not allow us to capture such effects, even 
at the free-action level. Indeed, had we kept them in the 
equations of motion, the dissipation terms we obtained 
in pop would have given rise to two friction terms (one 
per field) with sign opposite to the usual Hubble friction 
term. This anti-dissipation is therefore of a very different 



nature, and, just like the noise terms, they can be in- 
terpreted as coming from the constantly incoming modes 
through the window function into the coarse-grained the- 
ory. The fact that usual projection methods from the 
equations of motions miss this effect interestingly seems 
to suggest that dissipation arising from the fluctuation- 
dissipation theorem is a quantum effect (even though in 
the current case it turns out to be a subdominant one) 
which semi-classical methods, at least in their usual form, 
cannot capture. 

As a flnal note, we would like to emphasize that 
throughout this paper, mainly for the sake of clarity, we 
have been working with fluctuating fields on a fixed clas- 
sical background metric. In order to be fully consistent, 
one should explicitly calculate every noise term defined 
in terms of the slow-roll parameters, and keep gravita- 
tional corrections consistently up to whatever order we 
obtain. A priori, this might require expanding the equa- 
tions of motion for the mode functions up to second order 
in metric perturbations, which would make things more 
complicated. 



V. CONCLUSION 

In this paper, we have presented a heuristic but pre- 
cise derivation of the Langevin equations of motion for 
stochastic inflation. We then have presented an alter- 
native, more rigorous derivation from first principles of 
the formalism of stochastic inflation and pointed out the 
link between the presence of noise and the processes of 
decoherence and of classicalization. We then developed a 
self-consistent recursive method for obtaining a solution 
at consistent order in fields perturbations and in slow- 
roll, motivated by the case-study of multi-field hybrid 
inflation which we explore in a follow-up paper, where it 
has been suspected for many years that quantum effects 
and back-reaction can have a significant effects. 

In line with this, we developed a better understand- 
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ing of the loop expansion both in the bath fields and 
in the system-bath coupling within the CTP framework. 
In particular, we found how the term representing mode 
coupling found in the heuristic derivation emerges from 
the perturbative expansion. 

In the second paper of this series [3 If . we will apply 
our new recursive method to the specific example of two- 
field hybrid inflation potentials, as an example of how 
this method can be implemented in a generic multi-field 
context. We calculate the modified predictions emerg- 
ing from a consistently-implemented non-perturbative 
method for cosmological observables such as the CMB 
power spectrum. Most interestingly, we identifiy regimes 
of hybrid inflation where stochastic effects dominate over 
regular perturbative corrections. This result has poten- 
tially important implications for accurately constraining 
such small-field models of inflation. 

In future work, we wish to apply our method includ- 
ing mode coupling effects during the waterfall of hybrid- 
type models. This is of particular interest, since stochas- 
tic effects have been shown to completely dominate the 
classical evolution of the background fields trajectory in 



this phase [15|. Furthermore, there is plenty of scope for 
applying the method developed in the present paper to 
phenomenological models of inflation, an avenue which 
we plan to explore in future work. 
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APPENDIX A 



Path integral over the bath degrees of freedom : 

Leading order 

Using the bath-system splitting defined in sections |ll] and IIIIl as well as the Fourier modes decomposition of the 
bath degrees of freedom defined in p5|) , in the following our goal will be to perform an integration of the action over 
the bath degrees of freedom (the short wavelength fields) , 0> and tp^ , assuming the full fields $ and ^ are massive 
but have no interacting potential, in order to obtain an effective action 5off for the system of long-wavelength modes, 
and finally take the classical limit. 

In other words, we want to look at how to perform the path integral for the leading terms including bath fields, 
that is, the bilinear piece of the action. As we will see, this term will give rise to the noise terms we encountered in 
(fT8|) and p9|) . Using the notation we introduced in ((33|) - ([35]) of section IIIH we want to calculate: 



F[^,X] 



^yjdlpyj / V(l)y I X'^±e^^'^ba«'"'^bathJ 



^+ d7/.+ / p0±ei/'i'-[^^5A.0>+^^A, 



^>f^r>f 



VtP^e 



±pt / d*:l[5'A> A^V;>+X^A*'A>] 



(Al) 

where we are still using the matrix notation introduced in (j34p . Since A^ and A^ arc symmetrical operators with 
positive eigenvalues, we can then complete the square in the exponential and perform the Gaussian path integral over 
the bath fields: 



F['f,x] = Nexp 



1 i 

2h 



d^xd^x' ifi^ ix)A^iA^)-' A^^ix') + / dWx^(y)A^(A^)"'A^x(y') 



(A 2) 



where A^ is a normalization factor in which we absorbed factors of square root of determinant of A, and the arrows 



over A^ and A^ indicate the direction in which the time derivatives they contain are acting. Also, here, (A^) 
(A.^)~^ are Green's functions on the contour for the bath fields, (/)> and ?/>> respectively, in the CTP method: 



and 



ih{K^y 



:r0+(.T)0+(.T') 



0>(x')0j(x) 



m^ix)4>^ix')] 



ih{K^y 






i'>{y')i'tiy): 



r0;(y)0;(2/') 



(A3) 
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The entries of these propagators can easily be interpreted intuitively. The upper-left entry in each propagator repre- 
sents the amplitude for propagation between point x and x' on the "-|-"-part (or forward part) of the contour. The 
lower-right entry represent the amplitude for propagation if both x and x' are on the "— " or backward part of the 
contour. The off-diagonal entries represent the amplitude for a particle to propagate from an "cc" on the forward part 
of the contour to and "a;"' located on the backward part, and reversely. 

Writing the fields in Fourier space, and recalling that the bath fields are defined as in ((25|) . we can pull out of the 
propagator one factor of the window function on each side of each propagator: 



F['f,x] = Nexp 



2h 



dtdt' 






(27r)3 (27r)3 

T0u(t)0q(t')) (0k(t')0q(i) 
^' ^ /0k(t)0q(i')) {f4>v{t)4>^{t' 

X exp [<^ o X, ^ ^] , 



eaH 



Wh (^) l^{t')^-^{t') 



(A 4) 



where in the last line, we mean repeat the previous exponential term, but with all ip and (p fields replaced by x and f/; 
fields. Now, recall that the system fields <f and x are defined to be non-zero precisely on the long wavelengths where 
the window function Wh was chosen to vanish. The orthogonality (in the W ^ 6 limit) of ipk and Xfe with Wh means 
that the only non-zero terms in the above will be the ones where at least one time derivative is acting on all Wh^s. 



N exp 



2 \h 



dtdt' 



Nexp 



■ ■ (2^(2^" W^-W^*'^^ 
a^{t)^l^{t)Pt X 



'Mt)Mt')) (f0k(t)0q(t'); 



X exp [<y9 o X, f/* "^ V'] 

'' '• '• d^k 



2 \h 



dtdt' 



(27r)3 

oit - t')MtHUt') + 9{t' - t)Mt')K{t) 



N exp 



Mt)Kit') 

X exp [</j -r-> Xi f/* "H> V] 
2 U, 



e{t' - t)Mt)rw{t') + o{t - t')Mt'Mit) 



d'^xd'^x'a^t)^' {x)Pt 



(2^ 



0{t - t')4>i,{t)e-'^-'^(t>l{t')e'^-'^' + 9{t' - t)4>^{t')e'^-''' (l)l{t)er^^-'^ 

-0k(i')e'''-">k(i)e"''"" 
0{t' - t)cj)i,{t)e-'^-''<j)l[t')e'^-''' + e{t - t')Mt')e'^-''' (j>l{t)e-'^-- 



X exp [(f -n- X, <~> V'] 



a3Pt'a''{t')ip.qit') 
(A 5) 



Pt'a'(i')<^k(t') 



(A 6) 



Pt'a^it'Mx') 



(A 7) 



where from (|A 5P to (|A 6p we have used the operator Pt defined in psp as well as the expansion of the quantum fields 
in terms of creation and annihilation operators in order to evaluate the expectation values, 0k = ^k^k + 0k^k and 
similarly for ip\^. Also from ()A 6|) to (jA 7} . we have re-expanded the fields (f back in real space. 

The explicit expressions for the propagators above demonstrate that the upper and lower contours are not indepen- 
dent. In fact, the Green's functions written above contain a certain degree of redundancy. For example, we have that 
the trace of the matrix above equals the sum of the off-diagonal terms. Therefore, defining the linear transformation: 



U 



_ ( \/2 1/2 



1 



-1 



(A 



we obtain an equivalent representation of the bloc propagator we were working with. In this new basis, called the 
Keldysh representation, the system fields take the form we used in ([5^ to define the quantum and classical fields: 



(P =U'^= y 



- I ^c 



X -X 



- I Xc 
Xq 



(A 9) 
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where tfci Xc and Lpq, Xg define the classical and quantum fields, respectively [17j . In terms of these new variables, and 
making use of the identities Re[a&*] = {ab*+a*b)/2 and Im[a6*] = {ab* — a*b)/2 we obtain ([55)) after a short calculation. 
We therefore find the real and imaginary parts of II^ .^(x, x') contributing to the influence functional, which we can 
interpret as dissipation and noise kernels coming from the bath degrees of freedom, respectively, as discussed in section 
mil However, in order to be able to properly interpret the imaginary part of the influence functional as noise, we now 
proceed to a few aesthetic manipulations on the operator ([37]). 

First we rewrite Pt in a way that it contains only first time derivatives of the window function Wh'- 



a^it)[PtMt)]e-''''' 


-{0. 


aHt)WH 


+ 2a''{t)WHdt}Mt)e'''''' 




= dt 


a\t)WHMt)e-'''''' 


+ aHt)WHdtMt)e~'''-'' 




= dt ' 


a\t)WHMt)e~'''''' 


+ a^{t)WHlpM + q,p 



k 



a{t)H{t) 



{t)e 



nu[t)e 



(A 10) 



In the last line, for the mode functions of the quantum field we have assumed the form (j)k = {Pipit) + q^ [k/{aH)]} (f>]^ 
with some function of time q which contains all the fc-dependence through the combination k/aH = — fcr(l — e) (to 
first order in slow-roll, where r is, as before, the conformal time, and e = —H/H" is the first slow-roll parameter) 
and some function of time p with no fc-dependence. This is the case whenever mode functions can be written under 
the form of a Hankel function of the first kind, which is what we will be assuming in the following. For concreteness, 
in the case of single-field infiation, for a field with mass m where we neglect metric perturbations, we find: 



p{i)+q{-kT) = -H{t) 







(All) 



where v^ assumes the standard value of | — ^ -|-3e (to first order in slow-roll). In general, the exact functions p and q 
will be different for the two fields </> and ^ we are considering here, and will also depend on the level of approximation 
we are making for solving the equations of motion of the mode functions. 

Once substituted back in ([36)) . the first term of (|A lOp can be integrated by parts to yield (recall that, as explained 
in the main body of the text, the second term of (|A 10[) gives rise to both a mass-renormalization term and dissipation, 
which are negligible in the context of slow-roll infiation [50, l53] ) : 



SlK = 



d'^x<rx'a^{t)a!^{t') 



d^k 



2hJ "■ '"'" '" 'J (27r)3 

X { [Plit') + q; i-kr')] v,{x') ~ drv,{x')] 
+ (X ^ <^) , 



{[p^{t) + q^ (-fcr)] ^,{x) - dtVg{x)} Re I W{t)<j>^{t)W{t')cj>*^{t')e-'''<^-^"^ 

(A 12) 



as argued in section [Til Wh is constrained to have the form Wr {j^)- Replacing Wh = 

~k/{ae) 1 1 + H / H^ J {dk./(eaH)WH) ~ Hkr (1 — e) dkrWn (where the last equality holds to first order in slow-roll), 

we can perform the angular part of the integral over cPh and rewrite the above expression more conveniently in matrix 
form: 



SiA = ^ d'xd^x'a^t)a^t')[a^{x),^,{x)]jV^ix,x') 



at{x') 

iPq{x') 



+ (X ^ V') , 



(A 13) 



where we have defined (t'1^{x) = [p4,{i) — dt] fqix) and where the operators A^-'^,(a;,a;') are defined in ([4T|) . Since the 

j and t — > t' , we obtain 



two matrices M^"',,^(fcT, fcr') defined in (I42[) are Hermitian under simultaneous exchange of i 



<f,,ll!\ 



that the two Re 



]V[V,(fcr, fcr') are symmetric, and so are the two AV, (x,x') operators. 



This term gives an imaginary part to the effective action, and therefore cannot be properly interpreted as an ordinary 

-ax^/2 



part of the action. Using the Hubbard-Stratonovich transformation [67|,[68[, e 



VW^)IZ e-y/^^^^-"ydy, 



we introduce an auxiliary field y obeying Gaussian a-statistics and regularly coupling to the initial variable x. One 
can generalize such a formula to symmetric operators and introduce two real classical random fields per field in the 
system ^f , ^2 a-nd ^{, ^2 > each obeying the Gaussian probability distribution ([^0]) . and thus rewrite the influence 
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« c(l) 



functional as follows (neglecting an overall normalization factor): 
F[(p,x] =exp 






EE I ve.VitV [ef,ell exp ^ W dxa^t) [<j*{x), ^,{x) 



Cf(x) 



exp<{ - / d.Ta=^(t) [cr^(x), ^,(.t)] 






- / dxa^{t) [a'^^ix), ^,(.t)] 



I?er^e2''exp[0O7^] 



exp<{ - / dxa3(t) [cr^(a;), X9(a;)] 



ef(x) 



ef(x) 



-0 ^0 -i/; -i/; 
^1 J^2 'SI 5^2 



(A 14) 
(A 15) 



and the corresponding average ( ) ^ * ^ * c^i c^ ■ 

Sl ^S2 >?1 >?2 



where the two last equations define the probability density V Ci 7^2 

We therefore see how the influence action can be understood as adding an extra contribution averaged over classical 

noise sources to the total effective action. 

We are now in a position where we can go back to the total effective action. Writing it in terms of the quantum 
and classical fields instead of the fields on the forward and backward parts of the time contour using ([55]) 



V rx 



* 0(1) 



Vip^ 



Vx cxp 



- 5+„ - 5. 



(1) 



^ \^'-'sys "-"sys ^ "^lA 



Vip^-'Vx"^' I V^,V^2V [$1,6] exp ( ^ / d^xa^it) {ip^ (n - ml) ip, + x, (□ - O Xc 



-vWc + '^.Xc + 



Xq 



).V[ 



"Pq Xq 



[cr^(a;), ifiqix)] 






+ KW, Xqix)] 






(A 16) 



In order to explore the classical limit of the effective action, we proceed to the field rescaling fq,Xq ~^ f^^qi ^Xq- 
With this redefinition, a first order expansion of the effective action in tupq, hxq corresponds to the ^--independent, or 
classical, sector of the action, which allows us to obtain (|43p . 



APPENDIX B 



Perturbative Expansion 



In this appendix, we define the perturbative expansion in terms of the usual sum over loop diagrams of increasing 
order, which allows us is section IV B-C to perform the path integral over the bath degrees of freedom up to one loop. 
Using the notation we introduced in ([55)) - ([55)) of section IIVI we can define the free bath action on the forward and 
backward time contours: 



'hi bath bathj 



,i/d^^i[05A<^)0>+V'?:A<^'V'>] 



(Bl) 



In order to define this perturbative expansion [20|, |6l|, we introduce four currents, one for each field on the forward 
part of the time contour, and one per field on the backward part. Using the vector notation Jj ^ — [Jt.,/,, ~ ^<i^i/,]j 
we obtain the modified, unperturbed influence functional: 



F^'^J^J.^.] 



d(t>t,d^t, 



2?0± / VtPl 



xe 



i K+,-sL-t.+/.? d^-Jii>+!:! d^-Jii'>\ 



(B2) 
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Defining the average of an operator depending on the bath degrees of freedom Q[4>y, V'>] ~ Q [(/)>, V'>> 0>7 '0>] • 



Q[<f>>,i': 



d(l)+^di'+^ / V(t>^ / X'V>e^^ -='"■" ''-"'JQ[(/>>,i/' 



Q 



h 5 



tSJ+ixY iSJ+ix)' iSJ^ixY iSJ^ix) 

we can write the fuU influence functional F in terms of F^^': 
F[tf,x] = /e--fe[^bath(0t^J)-l4.th(0;,V'>)+V'.„t(v+,x+,0t'Aj)-vi„t(v-,x-,0;,'A>)] 



F(i)[j0,J^ 



J^ — J-ip—0 



cxp ■ 



y P S ^ S \ ~ 

V'bath — ^ -,1, , , — ^ ^1, , — Vbath 



ft (5 





h 5 



(B3) 
(B4) 

(B5) 



+V,,,t\^^,X^,- 



h 5 h 5 " 



I SJ^ (s) i SJ. (s) 



.y - - fi S h S 

-Vint [-^ ,X , T- 



zW+(5)'zW+(,s); \^'- '- ' ^SJ- ^5J^{x) 

Expanding the exponential, we obtain the perturbative expansion we are after: 

FW,X] = 
(l)n- 



F(i)[J^,j0] 



J4, — J^—0 



(B6) 






+ 



(v^2^(^+,X+>^,^^))^ + (v^?,t(¥'",X">;>>))^-2(vi„t(^+,X+,'/'^>^)V^nt(^-,X">;>>))J + 



where the ellipses represent terms with at least three powers of V . The terms already written down will give rise to 
terms similar to the ones we wrote down in (jlSp and (J19p and which we argued represented interactions between the 
quantum fields and the coarse-grained fields (that is, the terms that are on the fourth line of each equation). 



APPENDIX C 



Perturbative Calculation of the InDet Terms 



Our goal in this appendix will be to expand the determinants appearing in the expression for the free infiuence 
functional (j56p in a way that they can be explicitly calculated perturbatively. Doing so we will link this expansion 
to a diagrammatic expansion, in order to show that including the quadratic terms as part of the free propagators of 

the bath fields we get the same result as if we had treated the ^(v3='=)^(-0='=)^, ^ ix'^)'^ (4''^)'^ ^ and |(x''^)^(V'^)^ as 
perturbations and obtained their one-loop contribution to the perturbed infiuence functional. 

In order to evaluate the logarithm of the determinant of the operators defined in (j56p . work with the conformal 
time T in coordinates where the Fourier transform of □ = ^ [^^t + ^^] is the four-norm k^ /a^ In order to simplify 
the notation of the propagator bellow, we will set a = 1 for now, but to restore the factors of the scale factor in the 
calculation performed bellow, one only needs to add a factor of 1/a^ every time a four-norm k^ appears. 

We start by rewriting it as a trace, which we can evaluate as the sum of their eigenvalues: 



-Trln(A 



(2)^ 



d^k 

(2^ 



In [{^k' - ml - g^ix+f) {k^ + m| + g^x'T) - • O] 



(CI) 



Here the fi factor sitting at the front of the four-momenta integral represent the four-dimensional volume over which 
we are evaluating the path integral. It arises in this form since for simplicity we have written the x fields as if they 
had no time or spatial dependence, which is the case over the k-scales we are considering in the momentum integral, 
but once we re-introduce their spacetime dependence, this factor will transform into a four-dimensional spacctime 
integral. 
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Also, we have kept the zero contributions from the off-diagonal terms of A^ ' because wc wish to re-write this 
expression in terms of bath propagators. Specifically, the Feynman (Dyson) propagators lA^j^j^i) and Wightman 
propagators ^A^ > satisfy the inhomogeneous and homogeneous Klein-Gordon equations, respectively: 



□ . 



g'ix 



+(-)^2 



i^<t,,F(D){x,y) = {-)i5'^{x - y) 



□ , 



g'ix^^-^ 



iA^^<(»(a:,y)=0.(C2) 



We therefore replace the two '0' in the above equation ()C 1|) by the appropriate product of the two Wightman 
propagators and their kernels, recalling that they come form the explicit calculation of the determinant. We want to 
expand this expression around the constant free, massive Klein-Gordon propagator (free of system fields dependence), 
therefore we factor out those constant leading factors and collect them as follow: 



1 



Trln(A 



(2)^ 



d^k 



In 



(2^)4 

(-fc^ - m|)(fc^ +m|) 



>)(fc^ 



-fc2 



-fc2 - 



-g^x 



-H\2 



fc2 



-g'{x-f 



-g [x 



iA0,<(fc) 



-g [X 



iA0,>(fc) 

(C3) 



Looking at the Wightman propagators appearing in the above, we would like to write them as an expansion 
around the free, massive, system fields-independent ones, A\ l.(fc). Since A^>{k) are each proportional to a delta 
function in Fourier space, and therefore their Fourier expansion should be interpreted as a requirement to integrate 
by parts the Fourier integral appearing at the front of the above equation, once we perform that expansion we obtain 
^ri,.p{k) = A\ >(fc) -I- O^g"^) plus other terms which will give zero as a c-number once we have imposed the fact that 

free Wightman propagators are solutions to the free homogeneous Klein-Gordon equation: (fc^ -f m'^)iA, >(A:) 
Simplifying, we obtain: 
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(fc2- 

(2^l(-fc2-m|) (fc2+m|) (-fc2-m|) (fc2+™|) 
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where in the last line we have used our assumption of the validity of perturbation theory around small coupling 

constants to Taylor expand the In's up to 0{g'^) and where A\ [as defined in (|34l) and ([35]) ] is the integration kernel 
for a free, massive scalar field and so has no system field dependence, which means the terms we have collected in a 
new Tr In term boil down to being simply an overall constant multiplying the influence functional. 

Next we want to re-write this expression by introducing factors of the free bath propagators. To do so, we use the 
fact that free propagators are solutions to the Klein-Gordon equation: 
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(C4) 



We obtain 
1 



-Trln(A 



(2)^ 

4> > 
2 



1 



Tr ln(A 



(0)n 

> 



d^k 

(2^ 



.(0) 



.+ A2a(0) 



g' (AS:kfc)(x 

(0) ^^^^,,+^2A(o) 



^0,-D 



{k){x 



(C5) 
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+ . 



25 

As advertised at the beginning of the appendix, this result matches exactly the contributions coming from one-loop 
diagrams up to 0(1?^) (which is equivalent to O [(0^)^] in terms of the bath degrees of freedom) we would have 

obtained if we had treated the ^(0^)^(x^)^ term in the action perturbatively instead of as a time-dependent mass 
for the bath fields. Indeed, the first two terms in the above expression represent contributions from the following two 
diagrams: whereas the remaining terms are the contributions coming from the following 0{g ) (which is equivalent 



/I 



11 + 



X 



+ 



X 



+ 



II II 

-vv //- 

sLoj£^ 



X 



X 



to O [{(f)^)"^] in the bath degrees of freedom) diagrams: 




Reintroducing the spacetime dependence of the system fields, and rewriting the propagators in position space and 
in matrix form, we obtain: 
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-Trln(A;-^^) = --Trln(A^"') - -'- / d^x [ x+(^), X'ix) ] 
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(C6) 



Rotating to the Keldysh basis, we finally obtain (|72p . 

To derive (j73|) . one proceeds in an analogous manner, however one obtains more contributions at every order in 

the couplings, since there are two time-dependent mass terms in the non-zero entries of A\ instead of one. The 

additional diagrammatic contributions we obtain are, at 0{g^) or 0{\) in the couplings (i.e. O [(V'^)^] or O [((/>^)^] 
in the bath degrees of freedom): 







^ 



^ 



while the terms coming from the expansion of the InDet(AL') ^ term giving contributions at 0{g'^) or 0{\'^) (or 
equivalently O [(■0^)'*] or O \{k)'^'\) ^re equivalent to the contribution coming form the following diagrams: 
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